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ee temperature dependent material traits on random geometries. 


Available online: 24 Aug 2021 Further, the calculations can be pooled with structural analysis 
© 2021 The author(s). Published by AI Publication. plus convection can also be taken in report. Thermoelectric cooler 
This is an open access article under the CC BY employs Peltier effect for dissipating heat in an electronic casing 
license structure. It shows exceptional rewards over conservative cooling 
(https://creativecommons.org/licenses/by/4.0/) skill via quiet process, extended life span, and effortless integration. 


, ; Nevertheless, Joule heating results in the accumulation of internal 
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i ; ; heat thereby exposes thermoelectric cooler towards the risk o 
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, thermo-mechanical breakdown all through continuous operations in 
thermoelectric generator. , , , 

pragmatic thermal surroundings. In this paper, a 3D module of 
thermoelectric material MoSi2 is designed on the way to examine 
the thermoelectric effect of the material taking into consideration 
the temperature reliant TE material traits. The transient behavior is 
also observed. The results can be openly used intended for 


consistent design considerations and optimized thermoelectric 





devices in engineering. 


I. INTRODUCTION coefficient, o indicates the electric conductivity and À 


The thermoelectric effects within finite element analysis indicates the thermal conductivity. Generally the material 


(FEA) can be realized by means of the COMSOL- traits rely on the temperature moreover may perhaps be 


Multiphysics platform. It lets thermoelectric calculations anisotropic. At this juncture simply isotropic substance 


among temperature dependent material traits on random traits are worn. For anisotropic resources, the appropriate 


geometries [1]. The field equations in thermoelectric matrices are taken in consideration. The transient 


coupled intended for temperature as well as electric magnetic fields are also not taken in consideration. The 


potential under steady state calculations are described as projected equations are as a consequence to the coupled 


5 equations in [2] or the text referred within [3]. 
-V((oa?T + A)VT) — V(oaT WV) = o((WV) + 


aVTVV) (1) II. GEOMETRICAL MODEL 
and COMSOL Multi-physics allows the execution of ordinary 
V(oaVT) T V(oVV) = 0 (2) random partial differential equations (PDEs) intended for 


ae the field variable u over a one to 3D section Q. Two PDE 
where the material traits a indicate the seebeck- B . = 
modes are worn: The “Coefficient-Form” as well as the 
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“General Form”. 


2 
Cad 4 (—=cVu — au + y) + b. Vu + au = f 


a at 
(3) 
The thermoelectric field equations at this instant are 


altered into the “coefficient form” as follows. In the midst, 
the vector value of the field variable is defined by 


s (T 
u=(_,) (4) 
the coefficient c in (3) is 
A+ ca*T caT 
5 
( oa o ) ©) 


Intended for transient calculations the capacitive influence 
need to be neglected. Generally it is satisfactory to mull 
over merely the thermal capacity (heat capacity C, density 
p). Then d in equation (3) is 


d = (") (6) 


The subsequent examples show outcomes of calculations 
for characteristic thermoelectric applications. The material 
traits for the calculations with temperature independent 
values are depicted in table 1. Here characteristic values 
for Molybdenum Silicide MoSiz were taken from [4] and 
Copper was taken from [2]. Temperature dependent 
material traits were interpolated by means of cubic splines 
(figure 1-3). 
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Fig.l: Temperature dependent Seebeck coefficient of 
MoSi2 and cubic spline interpolation. 
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Thermal conductivity (Vink) 
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Fig.2: Temperature dependent thermal conductivity of 
MoSi2. 


Electric conductivity (10° AA) 
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Fig.3: Temperature dependent electric conductivity of 
MoSi2. 


Table.1 Numerical material properties. [4] 


Material Properties 
Thermal Conductivity 66.2 W / (m.K) 


Electric conductivity 3.28e° S/m 
Seebeck Coefficient 3.9e° V/K 
Heat capacity at constant pressure 430 J / (kg.K) 
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Table 2: Temperature dependent material properties. 


5 
fren | 





Hl, THERMOELECTRIC COOLER 


The geometry of a straightforward thermoelectric cooler 
comprises of solo p-type semiconductor component with 


dimensions | x 1 x 6 mm3. It is sandwiched by two copper 
electrodes of 0.1 mm in thickness (Figure 4). 
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Fig.4: A p-type thermoelectric element. 


The base is kept back at temperature 300 K along with OV 
of voltage. At the top of the upper electrode, a current of 
0.7A was applied. The resultant distribution of 
temperature is revealed in the middle. A temperature 
difference of nearly 70 K is achieved. Table 1 shows the 
(constant) material properties. Figure 4 shows the result of 
the calculation. In the center, the temperature distribution 
shows that the cold side temperature is at 230K. The 
associated voltage is shown right. To drive the current, a 
voltage of 50 mV is needed. 


IV. TRANSIENT OPERATION 


Figure 5 shows the outcome of a time reliant computation. 
The chart reveals the transient cold side temperature with 
temperature dependent material parameters. The short 
current pulse leads to a momentary temperature plunge of 
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about 3K. Such super cooling effects are also described in 


[5]. 
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Fig.5: Transient calculation of Peltier super cooling. 
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A tiny current pulse leads to momentarily lesser 
temperatures at the cold end. In such transient 
computation, barely the thermal capacities as suggested 
by equation (6) in the midst of the heat capacities in 
addition to densities are represented in table 1. 


V. _THERMOELECTRIC GENERATION 


In order to simulate a thermoelectric generator, the earlier 
mentioned semiconductor component was worn yet again 
by means of the changeable material traits (figure 1 — 3). 
The top side of the higher electrode was adjusted to 373K, 
whereas the base of the lower electrode was adjusted to 
273K along with OV. Figure 6 displays the outcome of the 
current - voltage characteristics of the thermoelectric 
material and Figure 7 displays the outcome of the current 
power characteristics of the material. 
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Fig.6: Current-voltage characteristics of the 
thermoelectric material 
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Fig.7: Current- -power characteristics of the 
thermoelectric material 
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In accordance to the properties, it was observed that the 
open circuit voltage of the component is computed to be 
about 2lmV, whereas the short-circuit current is 
computed around 220mA. The highest power output is 
observed as 1.22mW. 


VI. SUMMARY 


The accomplishment of the thermoelectric field equations 
using COMSOL multi physics 5.2 is projected. 
Thermoelectric computations may perhaps be finished for 
arbitrary geometries too. Anisotropy (not revealed here) 
as well as temperature reliance of the materials can also 
be incorporated. In addition, transient computations were 
made. It is probable in adding the structural analysis or 
convection effortlessly (not exposed here). 
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